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ABSTRACT 

This study was an evaluation of a radiational scheme 
for inclusion in a numerical weather prediction model. 

This scheme uses empirical formulations for atmospheric 
absorptivities, scattering and reflectivity and cloud 
and earth surface reflectivity to compute the solar in- 
solation received and absorbed by the earth and key atmo- 
spheric layers. Longwave cooling by the earth is also 
treated with empirically derived emissivities for water- 
vapor and carbon dioxide. This radiation model includes 
a two-layer cloud parameterization scheme and calculates 
the effect of an additional cloud layer on previous 
studies . 



In testing 


the 


validity of 


this mod 


el, use was 


made 


of Fleet Numeri 


cal 


Weather Central data 


in the form 


of 


vertical soundi 


ngs 


on a spring 


day over 


the ocean. 


This 



data was used to calculate the radiation balance at the 
tropopause, the earth's surface and the key layers of 
the atmosphere. 



2 



TABLE OF CONTENTS 



I. INTRODUCTION ----------------- 10 

II. DATA PREPARATION --------------- 13 

A. TEMPERATURE TERMS- ------------ 13 

B. MOISTURE TREATMENT IN THE SOUNDING - - - - 16 

C. PRESSURE-SCALED ABSORBER MASSES- ----- 19 

D. CLOUD PARAMETERIZATION ---------- 21 

III. SOLAR RADIATION- --------------- 22 

A. PARTITION OF SOLAR INSOLATION- ------ 22 

B. DISPOSITION OF F(S) INSOLATION ------ 24 

1. Clear Sky Case ------------ 25 

2. Cloudy Sky Cases ----------- 26 

3. The Weighted F(S) Insolation ----- 27 

C. DISPOSITION OF F ( A ) INSOLATION ------ 27 

1. The Clear Sky Case ---------- 27 

2. Overcast in both Cloud Layers- - - - - 30 

3. Disposition of F (A) Insolation with 

an Upper Overcast- ---------- 33 

4. Disposition of F ( A ) Insolation with 

a Low Overcast ------------ 33 

5. Weighted F (A) Surface Insolation - - - 36 

6. Atmospheric Layer Absorption of F(A) - 37 

a. Overcast in both Cloud Layers- - - 37 

b. Upper Overcast Only- ------- 37 

c. Lower Overcast Only- ------- 38 

d. Composite Weight ed Absorption 

by Layers- ------------ 38 



3 



e. The Absorptivity (ABA) by Layers - 38 

D. REFLECTION OF F (A) AND F(S) ENERGY - - - - 39 

1 . Albedo ---------------- 39 

2. Composite Absorptivity ( ABG ) by Earth- 42 

3 . Check on Computations- -------- 42 

E • STATISTICAL COMPUTATIONS --------- 42 

1. Clear Sky Cases- ----------- 42 

2. Relationship between Albedo in the 

Cloudy and Clear Sky Cases ------ 44 

3. * Statistical Relationship between 

Ground Insolation in the Cloudy and 
Clear Sky Cases- ----------- 46 

4. Absorptivi t ies in Air- -------- 47 

5. General Conclusions- --------- 47 

IV. TERRESTRIAL RADIATION- ------------ 48 



A. 

B. 



C . 
D. 
E . 
F . 



THEORETICAL AND EMPIRICAL BASIS 
* 



NET FLUX F ^ Q , WITH CLOUDY SKIES ------ 

* 

1. Net Flux Clear Skies - - - - 

* 

2. Net Flux F 1 Qf with Overcast Clouds 

in both Layers - - - - - - - - - - - - 

* 

3. Net Flux F 10 f with Overcast Clouds 

in the Upper Layer Only- ------- 



4. 



* 

Net Flux F , with Overcast Clouds 
in the Lower Layer Only- - - - - - 



* 



5. Composite F^ Calculations ------ 

* 

NET FLUX F _ , WITH CLOUDY SKIES ------ 

6 

* 

NET FLUX F , WITH CLOUDY SKIES ------ 

APPLICATION TO HEAT BALANCE COMPUTATIONS - 
STATISTICAL RESULTS- ----------- 

1. Downward Flux with Clear Skies - - - - 



48 

51 

51 

53 

53 

56 

56 

58 

63 

68 

68 

68 



4 



2 



. Statistical Relationship between Cloudy 

Sky and Clear Sky Cases- ------- 70 

3. F* (CL) Relative to F* (0, 0)- - - - - 71 

V. MERIDIONAL CROSS-SECTIONAL DEPICTION OF 

RESULTS- ------------------- 74 

VI. CONCLUSIONS- ----------------- 82 

LIST OF REFERENCES ----------------- 84 

INITIAL DISTRIBUTION LIST- ------------- 87 

FORM DD 1473 -------------------- 88 



5 



LIST OF TABLES 



I. a. Form of Original FNWC Grid Soundings- - - - 15 

b. Form of the Corresponding Radiative 

Soundings ----------------- 15 

II. Total Extraterrestrial Insolation (FAD J ) at 
Each Grid Point of the Three Data Lines 

(ly min-1)- - - 40 



6 



LIST OF FIGURES 



1. Five-Layer Radiative Sounding Model Used 

in this Study ----------------- 

2. FNWC Polar Stereographic Grid and Meridians 
(1,2,3) Selected for Study- ---------- 

3. Schematic Representation of F (A) Flux 
Disposition in the Clear Sky Case ------- 

4. Schematic Representation of F (A) Flax 
Disposition in the Case of Two Overcast Layers- 

5. Schematic Representation of F (A) Flax Dis- 
position with an Upper Overcast Layer Only- - - 

6. Schematic Representation of F (A) Flax Dis- 
position with a Lower Overcast Layer Only - - - 

* 

7. Terrestrial Net Flux F 

a . Clear Sky Case- -------------- 

b. Two Overcast Layers ------------ 

c . High Level Overcast ------------ 

d . Low Level Overcast- ------------ 

* 

8. Terrestrial Net Flux F 

6 

a. Clear Sky Case- -------------- 

b. Two Overcast Layers ------------ 

c. High Level Overcast ------------ 

d . Low Level Overcast- ------------ 

* 

9. Terrestrial Net Flux F 

a. Clear Sky Case- - - - - - - - - - - - - - - 

b. Two Overcast Layers ------------ 

c . High Level Overcast ------------ 

d . Low Level Overcast- ------------ 



11 

14 

28 

31 

34 

35 

52 

54 

55 

57 

59 

60 

61 

62 

64 

65 

66 

67 



7 



79 



10, Key to Mean Meridional Cross Sections 

11. Mean Meridional Cross Sections 

a. Tropical Section ------------- 80 

b. Higher Latitude Section- --------- 81 



8 



ACKNOWLEDGEMENTS 



The author wishes to express his appreciation to 
his wife for her support and for typing the rough draft 
of this thesis. 

Appreciation is also expressed to the author 1 s 
thesis advisor, Professor F. L. Martin, for his sugges- 
tions, advice, guidance and support in this research. 



9 



I 



INTRODUCTION 



This study was an attempt to produce an empirical 
radiative parameterization for inclusion in the Fleet 
Numerical Weather Central (FNWC) primitive equation pre- 
diction model. 

The radiative parameters described here include two 
cloud layers of fractional amounts CL ( 1 ) and CL(2) at 
each gridpoint instead of a single cloud layer as pre- 
viously used (Plante, 1973). The possibility of an addi- 
tional cloud layer tends to add substantial complexity to 
the radiation physics, but in general this paper repre- 
sents an extension of that described previously by Martin 
(1972). Some of the details of this radiation model are 
also based upon an unpublished manuscript (Martin, 1974), 
which in turn has some similarities to the multi-cloud 
models of Arakawa (1972) and Rodgers (1967). 

The vertical scale used in this study was based on 
the five sigma (a = p/TT) levels of the FNWC primitive 
equation model. The levels in this study are called k- 
levels (k = 1, e .., 10). The k-levels are used for the 

purpose of data input as shown in Figure 1. Cloud cover 
is dealt with by means of a two-layer parameterization, 
with one high-level cloud layer from k = 4 to k = 6 and 
a low-level cloud deck from k = 8 to k = 9, in amounts 
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Figure 1. Five-layer radiative sounding used in this 
study. Levels are identified by their values on the k- 
scale, while layers are identified by their level boundary 
indices (e.g., (8, 10) in parentheses. Pressure-scaled 

water vapor and C0 0 mass increments AU and AC, respective- 
ly, as well as mixing ratio (q) are introduced at the odd 
levels while temperature is introduced at all levels. 
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CL ( 1 ) and CL(2), respectively. In the computations for 
pressure-scaled water vapor and CC^ absorber masses a 
realistic method was used to extrapolate water vapor mix- 
ing ratio values above the level of k = 4, which is nor- 
mally the highest moisture level available from radiosonde 
data. The system for computation of the net flux of ter- 
restrial radiation makes use of new empirical absorptivi- 
ties of water vapor and CO (Sasamori, 1968), which are 
incorporated, in the Martin (1974) manuscript. 

To evaluate the radiative heat transfers of this study, 
three meridians derived from FNWC analyses were used. 

These meridians are all from the Pacific Ocean area and 
most of the radiative soundings were in the Northern Hemi- 
sphere . 
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II 



DATA PREPARATION 



The data used in this study was in the form of three 
meridians of 25 vertical soundings over the Pacific Ocean. 
These meridians were located at 125°W, 170°W and 145°E as 

shown in Figure 2. 

The soundings on these data lines were taken from the 
original FNWC analyses of temperature at the surface and 
at nine other standard levels up to 100 millibars for 
0000GMT, 25 April 1973. The six moisture values were given 

by the vapor pressure at sea level, and by dew-point de- 
pression for levels from 925 mb to 400 mb. Table 1(a) 
gives as an example the original sounding at gridpoint 
(1, 1) in the FNWC map. The data-proc e s s ing methods 

applied to the temperature and mixing ratio data for the 
determination of the water vapor and CC^ absorber masses, 
and the cloud amounts within required k-level boundaries 
(Figure 1) are detailed in the following subsections. It 
should be noted that all soundings start at sea level, 
and as a reasonable approximation may be considered to 
have surface pressure 7T = 1000 mb. Hence the k-levels of 
Figure 1 became 1000, 900, 800, 700, 600, 500, 400, 300, 

200, and 100 mb, respectively. 

A. TEMPERATURE TERMS 

The temperatures were listed at each mandatory level 
between 1000, ..., 100 mb. The radiative sounding 
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Figure 2. FNWC polar stereographic grid and 
(1,2,3) selected for study. The longitudes 
angles of the sun relative to local noon are 
each meridian. 



meridians 
A and hour 
shown for 
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TABLE 1(a). Form of original FNWC grid soundings. The 
humidity parameter between 925, ..., 400 mb is dewpom 

depression . 



Pressure (mb) T (°C) 

1000 28.5 

925 20.0 

850 16.2 

700 7 ‘ 1 

500 -9.6 

400 -21.4 

300 -36.4 

250 -45.3 

200 -55.5 

150 -66.9 

100 -80.3 



Humidity Parameters 

27.7 mb 
3 . 2 ° C 
2 . 9 ° C 

5 . 0 ° C 
4 . 7 °C 

2 . 0 ° C 



TABLE 1(b). Form of the corresponding radiative sounding 
with temperature and mixing ratio listed at k-levels. 



Pressure (mb) T (°K) 

1000 299.0 

900 291.7 

800 286.6 

700 280.3 

600 272.6 

500 263.6 

400 251.8 

300 236.8 

200 217.7 

100 192.9 

0 192.9 



Mixing Ratio 

18.35 g(kg) 
16.50 
12 . 97 
7.84 
4.75 
2 . 97 
2.21 
0.73 
0.23 
0.04 
0.00 
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temperatures at certain k-levels were obtained by inter- 
polation. The temperature T^ Q was set equal to the listed 
air temperature. At the top of the atmosphere, the tem- 
perature was assumed to be isothermal from 100 mb to p=0, 
with the mixing ratio also approaching zero at p=0. 

The interpolation scheme used was the three-point 
Lagrangian formula of the form: 

T k ■ T O <<I 'k- P l HI 'k- P 2 ),/<(P 0- P l )<I '0- P 2 ,) 

+ V<P k -P 0 > <p k- p 2 ) > /((P rV IP 1- P 2>> <2 - 1) 

+ V’VV (P k -P 1 ))/<<P 2 -Po ,<P 2- P l ) > 

Here lies above P^., where k is the pressure level de- 

sired for T -interpolation, and P ,P lie below level k. 

K X ^ 

For example, with k = 9, P^ and P^ were taken as 925 mb 
and 1000 mb, respectively. The parameters T^, T^ and T 
are temperatures at the same levels as P^, P^, and P 

respectively . 

B. MOISTURE TREATMENT IN THE SOUNDING 

The moisture parameters from the original sounding 
(Table 1(a)) were converted into the mixing ratios at each 
of the original sounding levels. At k = 10 the surface 
vapor pressure, e was used to calculate the mixing 

ratio by use of Equation (2-2) : 

= .62197 e . /1000 (2-2) 

^1000 air 

where e is in mb. 

air 
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In calculating the mixing ratios at the remaining 
original data levels Equation (2-3) was used (after p. 63, 
Fleagle and Businger, 1963) t 



q (P) = 



(A- — ) 
.62197 e _ 



exp 



-B (T-T ) 
t T T ~~ ] 



1000 



(2-3) 



Here q(p) is the mixing ratio at level p. 

A — 21.656 
B = 5418. 0°K 

T = temperature at level p (°K) 

= dew point at level p (°K) 

T-Tj^ = dew point depression (°K) 

The q-values were then interpolated to the k-levels using 



an equation similar to Equation (2 1) . The results of 
this interpolation are shown in an example. Table 1(b). 

Accurate humidity information is not generally avail 
able at or above 300 mb. To find the q-values at k = 3, 
2 and 1, an extrapolation formulation was used 



q.(£j = (_E_) 

' 500 ; 



(2-4) 



This formula is similar to that developed by Smith (1966) . 

However, instead of using the surface values of pressure 

and moisture in Equa tion ( 2- 4 ) use was made of the pressure 

and moisture values at the 500 mb level. Equation (2-4) 

was solved for the parameter X using the "best-fit" between 

the variables y = log 3— and x = log . This was 

q 5 
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done by using the six values of q of each original sound- 
ings corresponding to p = 1000/ 925/ 850/ 700 f 500/ and 
400 of Table 1(a). The "best-fit" values of X were then 
used in Equation (2-4) to extrapolate the q-values to 
p = 300, 200 and 100 mb. 

The usual formula for a best-fit parameter X relating 
y. and x.^ samples is given in standard texts on statistics 
(e.g.. Crow et al, 1960) 



X = 



i = 6 

E ' , i x i 

i = l 

\L 6 2 

E x i 

i = l 



A further statistic computed in connection with each X 
profile was the correlation coefficient 




The correlation coefficient for each X-profile was general- 
ly in the range .95 to .99, which tended to enhance the 
credibility of the usefulness of the X-profile technique. 

The mean values, X, together with the standard devia- 
tion, a (X) , based upon the 25 soundings on each meridian 
are listed below: 

Meridian 1 Meridian 2 Meridian 3 

X 3.143 3.712 2.931 

0 (X) 0.799 0.915 1.094 
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C. PRESSURE-SCALED ABSORBER MASSES 

With the mixing ratios computed at each odd k-level 
(Figure 1), it is possible to compute the pressure-scaled 
water vapor absorber masses in each layer. In Figure 1, 
these water vapor layer-masses were denoted by Au f but 
here we develop notation which is more representative of 
the layers under consideration* For example, the pressure- 
scaled water vapor mass in layer (8, 10) is given by 



^9 A P P 9 

U ( 8 , 10) = 



0.72 



g 1013.25' 



(2-5) 



where the pressure ratio raised to the exponent 0.72 (after 
Moller and Raschke, 1964) has been employed. For the layer 
(6, 8), an analog to Equation (2-5) was employed with q^ 

replaced by q^ and p^ by p^ . A similar procedure was fol- 
lowed in the topmost layers (0, 2), ( 2 , 4) and (4, 6). 
Equation (2-6) affords an example for the topmost layer 
( 0 , 2 ) . 



q, Ap 

u (0, 2) = [ 



1013 .25 



0 .72 



(2-6) 



In each of the five layers Ap is 200 mb. 

The algorithm for computing the integrated water vapor 
mass above the earth's surface is then as follows: 



u (10) 
u (8) 


= 0.0 
= U(10) 


+ 


U (8, 


10) 


U (6) 


= U ( 8 ) 


+ 


U (6, 


8) 


U (4) 


= U(6) 


+ 


U (4 , 


6) 


U (2) 


= U(4) 


+ 


U (2 , 


4) 


U (1) 


= U(2) 


+ 


. 5 (U (0, 2) ) 


U (0) 


= U(2) 


+ 


U (0, 


2) 
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The principle for the computation of the carbon dioxide 
scaled mass is similar to that for water vapor. Here, how- 
ever, it is customary to state the pressure-scaled mass in 

terms of the normal temperature and pressure (N.T.P.) 

2 

volume over a cm of the earth's surface. Thus 

C(8, 10) - 3.14 x < 10 lL 5 >°- 72 (2 - 8> 

s 

-4 

In (2-8), the fraction 3.14 x 10 is the constant mixing 

ratio of carbon dioxide by volume, while — — — is the N.T.P. 

^ ^s 

thickness of the layer AP = P - Pq. The quantity is 
the constant density in a standard homogeneous atmosphere 
for which T^ = 273.16°K everywhere and p^ = 1013.25 mb. 

With minor modification, Equation (2-8) becomes 

P 

C(8, 10) - 3.14 x 10- 4 (H) ( I3T f- IT ) 0 - 72 (2-9) 

RT 

H = — = 7.995 x 10 cm. 

9 

All other layer thicknesses of carbon dioxide are similar 
in form to C(8, 10), for example: 

0(6, 8) - 3.14 x 10- 4 (H)( 10 ^ -25 ) < 10 i 3 , 25 >°' 72(2 - 10) 

The integration for the layer masses of C(8), C(6), C(4), 

C(2), C(l) and C(0) follow exactly the procedure outlined 
for water vapor. 
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D. CLOUD PARAMETERIZATION 

In the calculations of the fractional cloud cover 
amounts in layers (4, 6) and (8, 9) the relative humidi- 

ties and thus the saturation vapor pressure at levels 
k = 5 and k = 9 must be available. The saturation vapor 
pressures at each level were calculated using Equations 
(2-11 (a)) and (2-11 (b)), and the relative humidity using 
Equations (2-12 (a)) and (2 — 12(b) ) s 



q^ (500) = 



q (900) = 

s 



(A- §-) 

. 6 2 197e 5 

500 

U- 2-> 

. 62197e 9 

900 



(2-11 (a) ) 



(2-11 (b) ) 



where A = 21.656 
B = 5418. OK 



RH(5) 



q (500) 
q s (500) 



RH (9 ) 



q (900) 
q s (900) 



(2-12 (a) ) 



(2-12 (b) ) 



The relative humidity values were then used to calculate 
the fractional cloud amounts at each layer using the fol- 
lowing formulation, which is similar to Smagorinsky 1 s (1960) 



CL ( 1 ) = 2 . 0 ( RH ( 5 ) ) - 0.7 (2-13 (a) ) 

CL ( 2 ) = 3.33 ( RH ( 9 ) ) - 2.0 (2-13(b)) 



21 



Ill 



SOLAR RADIATION 



A. PARTITION OF SOLAR INSOLATION 

In this study, following Joseph (1966), a solar con- 
stant of 2.00 ly min 1 was assumed. This is the value at 
the top of the atmosphere; however, it is assumed subject 
to a four percent attenuation by oxygen and ozone above 
the tropopause. Thus the effective solar constant assumed 
at level k = 2 in this study is 1.92 ly min ^ . 

The effective solar insolation at the tropopause was 
computed from 

F ( 2 ) = S — ] ~ 2 cos Z (3-1) 

R 



where 



Th 
in 
In 
va 
z e 



S = 
Cos Z = 
R/R 

m 

e ratio R/R i 
m 

the Smithsoni 

this study ba 

lue of R/R = 
m 

nith angle was 



effective solar c 

cosine of the zen 

ratio of the actu 
to the mean earth 

s listed as a func 
an Meteorological 
sed upon data for 
1.00606 was used, 
calculated by 



onst ant 

ith angle 

al earth-sun di 
-sun distance. 



tion of the 


Jul 


Tables (List 


, 1 


25 April at 


000 


The cosine 


of 



stance 

ian date 
958) . 
0GMT a 
the 



Cos Z 



Sin<{> Sin6 + Coscf) Cos6 Cos h 



(3-2) 



where 



(p = Latitude 

6 = Solar declination for April 25th = 12°56' 
h = hour angle 
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The value of the declination was taken from Table 169 
(List, 1958), and the value of Sincf) was calculated using 
two different formulas depending on the data-line used 
for the computations. These values of Sincf) along lines 
one, two, and three are calculated in terms of the FNWC 
map coordinates (I,J) as: 



Lines 1,3 Sincf) = 



Line 2 Sincf) = 



973.752 - 2 (32-1) 
973.752 + 2 (32-1) 

973.752 - (32-1) 

973.752 + (32-1) 



( 3- 3a ) 



( 3- 3b ) 



I = 1 , 25 

Here I is the abscissa distance on the FNWC grid. Figure 2, 
The soundings were all taken at OOOOGMT; therefore, the 
value for the hour angle h, can easily be calculated. For 
example, in this experiment solar noon existed at the 180th 
meridian and the hour angles h = 55°, 10° and 35° occurred 

along lines one, two, and three, respectively. 

Following Joseph (1966) a computationally simple parti- 
tion of the solar insolation was used for this study. 

Joseph divided the insolation at level 2, F(2), into two 

parts, one of which was subject to water vapor absorption 
but not to Rayleigh scattering. This portion is given by 

F (A) = 0.349F (2) ; (3-4) 
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the remaining portion of the insolation. 



F ( S ) = 0.651F (2) 



(3-5) 



is subject only to Rayleigh scattering (but not to water 
vapor absorption) in a clear moist atmosphere. The 
presence of two cloud decks in this study introduces 
cloud-reflectivity into both the F (A) and F(S) solar 
energy portions . 

The F(S) energy is composed of all wavelengths X < 

0.9 ym since here water-vapor absorption is negligible. 

The F (A) portion consists of those wavelengths X > 0.9 ym 
where absorption by water-vapor and carbon dioxide bands 
in the near infra-red is the dominant attentuation process 
in clear air. 

B. DISPOSITION OF F(S) INSOLATION 

Joseph (1966) found that empirical values of Rayleigh 
clear-sky albedo, (after Coulson, 1959), were best-fit 
using the least squares technique by the following formula 



occurs near 1000 mb. The value Sec Z is equal to (Cos Z) 
where Cos Z is calculated from (3-2). 

Another r ef lec t ive- type parameter which must be con- 
sidered in this study is the surface albedo a(G). 




(3-6) 



Here TT/p^ M in view of the fact that mean sea level 
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Over 



the land areas of the Northern Hemisphere a constant 
average surface albedo of 0.14 is utilized. However, the 
value over the oceans depends on cos Z. Since in this 
study all of the soundings were over the ocean, the 
equation 

a(G) = max {.06, .06 + .54 (.7 - cos Z) } (3-7) 

was used to calculate a(G), (after Gates et al 7 (1971)). 

Since a two-layer cloud model was used, four separate 
cloud cases are discussed below. 

1 . Clear Sky Case 

In this case the F(S) insolation is subject only 
to Rayleigh sky reflectivity, a(R), and the ear th- sur f a ce 
reflectivity, a(G), of Equations (3-6) and (3-7) 0 Con- 
sidering the possibility of multiple reflections between 
the earth and atmosphere, each of which diverts downward 
the fraction a(R) of the earth’s surface reflectance a(G), 
the clear sky insolation at the earth's surface after 
scattering is given by 

IS10(0,0) = F (S) [l-a(R) ] [l-a(G) ] / [1-a (R)a(G) ] (3-8) 

In the next three cases, where clouds were present, 
the cloud reflectivity values suggested by C. D. Rodgers 
(1967) were used. Thus for the higher level clouds (between 
k = 4 and 6), R(l)= 0.54 and for the lower level clouds (be- 
tween k = 8 and 9), R ( 2 ) is 0.66. 
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2. Cloudy Sky Cases 



We consider first the fraction of the grid area 
having clouds at both levels. The F(S) insolation ab- 
sorbed by the ground at such a grid point is given by 

IS10(1,1) = F(S) [(l-R(l)) (1-R(2)) (l-CL(G))] (3-9) 

* [1— ( (R(l)R(2)+R(2)0t(G)+R(l)0t(G) ) 

+ 2 (R(l)R(2)a(G) ) ] _1 

(after Arakawa, 1972) . The special symbols (1,1) denote 
CL(1) = CL(2) = 1.0/ which will temporarily be assumed. 

The other cloud-cover variations to be considered 
are those of CL(1) = 1.0, CL ( 2 ) = 0.0 and CL(1) = 0.0, 

CL(2) = 1.0, respectively. In the former case, R(2) = 

0.0 in (3-9) so that 

I S 1 0 ( 1 , 0 ) = F(S) (l-R(l) ) (l-a(G) )/(l-R(l)a(G) ) (3-10) 

In the latter case with clouds at the lower layer only it 
follows that 

IS 1 0 ( 0 , 1 ) = F(S) (1-R(2) ) (l-a(G) )/(l-R(2)CX(G) ) (3-11) 

Equations (3-8), (3-9), (3-10), and (3-11) were used to 

calculate the absorbed F(S) insolation at the earth's sur- 
face after only r ef lecti ve- type processes are considered 
to take place. However, in each case the sky was con- 
sidered to be either overcast or clear and an appropriate 
grid-area weighting scheme was necessary to correct to the 
actual cloud conditions CL ( 1 ) , CL ( 2 ) . 
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The Weighted F(S) Insolation 



The cloud combinations denoted by (0,0), (1,1), 

(1,0), and (0,1), respectively, actually occur with the 
areal coverage (or weights) at a grid point given by the 
following functions: 

W ( 0 , 0 ) = (l-CL(l) ) (1-CL(2) ) 

W ( 1 , 1 ) = ( CL ( 1 ) ) (CL ( 2 ) ) (3-12) 

W(1,0) = ( CL ( 1 ) ) ( 1-CL ( 2 ) ) 

W(0,1) = ( 1 - CL ( 1 ) ) (CL ( 2 ) ) 

at each grid point. 

When these weights were applied mu It iplicat ively 
in sequence with the right sides of Equations (3-8), (3-9), 

(3-10), and (3-11), then the composite surface-absorbed 
insolation, IS10 resulted. 

IS10 = IS10(0,0)W(0,0)+ IS10 (1 , 1) W (1 , 1) 

+ IS10 (1,0) W(1,0)+ IS10 (0, 1) W (0, 1) (3-13) 

Note that by subtraction of IS10 from F(S) we also 
arrive at that portion of the F(S) insolation reflected to 
space . 

C. DISPOSITION OF F ( A ) INSOLATION 

In the following sections, the fraction of the insola- 
tion which is absorbed by atmospheric water-vapor and 
carbon dioxide is discussed. 

1 . The Clear Sky Case 

In this connection the Manabe-Mo 1 ler absorptivity 

function 

0 3 0 3 

a(2,2k)= . 271 [U (2 , 2k) Sec Z] (3-14) 
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Figure 3. Schematic representation of F (A) insolation dis- 
position in the clear sky case. 
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was used. The form of Equation (3-14) indicates that the 
absorptivity a. is applied over the pressure-scaled water- 
vapor mass between levels 2 and 2k (see Figure 1) along 
the zenith slant-path angle Z. Therefore, assuming negli- 
gible water-vapor mass above level 2, the absorbed solar 
insolation in the layer (2,6) is 

A ( 2 , 6 ) = 0.271F(A) [U(2,6)Sec Z]° ,3 ° 3 (3-15) 

Similarly, the absorbed insolation in the layer (2,10) is 

A ( 2 , 10 ) = 0.271F{A) [U(2,10)Sec Z]° ,3 ° 3 (3-16) 

and the absorbed insolation in the layer (6,10) is given 
(see Figure 3) by 

A ( 6 , 1 0 ) = A ( 2 , 1 0 ) - A ( 2 , 6 ) (3-17) 

The direct transmission of F (A) insolation to the 
earth's surface was found by subtracting A(2,10) from F ( A ) . 
This direct transmission is partially reflected by the 
earth's surface (the value computed for a(G) in Equation 
(3-7) is applicable here) and the remainder, absorbed at 
the earth's surface, is given by 

IA10 (0 , 0 ) =F (A) {l-.271[U (2,10)Sec Z] * }(l-a(G)) 

(3-18) 

In the following subsections dealing with partially 
cloudy layers, it is necessary to adopt representative 
cloud reflectivities. For this purpose, values suggested 
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by C. D. Rodgers (1967) were used. These reflectivity 
values for the F ( A ) wavelengths differ from those listed 
by Rodgers for the F(S) wavelengths. We followed the 
procedure of considering overcast cloud conditions when- 
ever clouds occur and subsequently applying the weight- 
ing factors as was done in Section III(B). 

2 . Overcast in Both Cloud Layers 

Here, with clouds at both levels the cloud reflec- 
ivity values RA ( 1 ) = .46 and RA ( 2 ) = .50 were used. In 

addition Rodgers suggests cloud absorpt i vi t ie s of A(l) = 
.20 and A(2) = .30 (see Figure 4). With the following 

set of formulas, the downward directed insolations FA 4 ^ , 
etc., were calculated at each required level in order to 
arrive at a value for the insolation abosrbed by the 
atmospheric layers, and finally by the earth* s surface. 



A(2, 


4) ' 


= 


F(A) .271 [U(2,4)Sec Z]°* 3 ° 3 


( 3- 19a ) 


FA4 + 




= 


F(A) [ 1— . 27 1 ( U ( 2 / 4 ) Sec Z)° ,3 ° 3 ] 


( 3- 19b ) 


A(2, 


4) " 


= 


F (A) RA (1) . 271 [U (2,4) Sec Z]°* 3 ° 3 


(3-19c) 


A ( 2 , 


4) 


= 


A ( 2 , 4 ) ' + A ( 2 , 4 ) " 


(3-19d) 


A ( 4 , 


6) 


= 


FA4+A (1) 


(3-19e) 


FA6 + 




= 


FA4+ (1-RA (1) -A (1) ) 


( 3 - 1 9 f ) 
2 -1 


FA8 + 






FA6 + ( TD ( 6 , 8) ) [1-RA (1) RA (2) (TD (6 


, &)) ] 

( 3- 19g ) 


A (6, 


8) 


= 


FA6+ - FA8 i 


( 3-19h) 


A (8, 


9 ) 


= 


FA8+A ( 2 ) 


(3-191) 


FA9 + 




= 


FA8 + (1-RA ( 2) -A ( 2 ) ) 


(3-19 j ) 
2 - 1 


FA 1 0 1 




FA9 + TD ( 9 , 10) [1-a (G) RA(2) (TD(9, 


10) ) ] 

( 3- 19k ) 


A(9, 


10) 


= 


FA9+ - FA101 


cn 

o\ 
1— 1 
l 

00 
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Figure 4. Schematic representation of F (A) insolation dis- 
position in the case of two overcast layers. 
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where 



O 3 01 

TD(6, 8) = 1 - .271 [U (6, 8) 5/3] ° (3-20) 

TD ( 9 , 10)= 1 - .271 [U(9, 10)5/3] 0,303 

In the preceding equations A(2, 4), A(4, 6), 

A(6, 8), A(8,9), and A(9,10) represent the insolation ab- 

sorbed in the indicated layers. A superscript prime fol- 
lowing A(2,4) indicates absorption (by use of Equation 3-14) 
of the downward directed beam, while the double prime indi- 
cates absorption of the beam specularly reflected upward 
from the cloud top. The functions of Equation (3-20) repre- 
sent transmissivity after diffuse passage through a cloud 
or sub-cloud layer, which causes the effective replacement 
of Sec Z in Equation (3-14) by Sec Z = 5/3 (Katayama, 1966). 
The symbols FA6I, FA8I , FA9+ , and FAlO^ indicate the down- 
ward streams of insolation passing through the indicated 
levels. Note that the effect of multiple reflections be- 
tween clouds is incorporated in the denominator of FA8l, 
while that between the lower cloud and the ground has a 
corresponding term in the denominator of FAlCH, involving 
the reflectivities a(G) and RA ( 2 ) . The effect of the 
multiple reflectances as modeled here is to enhance the 
insolation at the lowermost reflecting boundary of the two 
under consideration. 

Finally, note that the, flux absorbed by the earth's 
surface under the conditions of an overcast layer at both 
levels is given by 

I A 10 ( 1 , 1 ) = FA10(l-a(G)) (3-21) 
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• Disposition of F(A) Insolation with an Upper 
Overcast 

In this case with only an upper overcast layer, 
the values of RA ( 1 ) and A(l) used are the same as those 
of the preceding section. However, a somewhat simpler 
set of disposition equations result: 

A ( 2 , 4 ) = F(A){.271[U(2, 4)Sec Z]°* 3 ° 3 

+ [1 —0271 (11 (2 , 4 ) Sec Z) ° ‘ 303 ] RA ( 1 ) } (3-22a) 

A ( 4 , 6 ) = F (A) A(l) [1-.271 (U (2 , 4)Sec Z)°' 3 ° 3 ] (3-22b) 

FA64-=F (A) [1- . 271 (U (2, 4)Sec z ) ° * 3 ° 3 ]' [ 1 -RA ( 1 ) - A ( 1 ) ] 

( 3- 2 2 c ) 

FA1 0 -V =FA6 + [1- . 271 (U (6, 10)5/3)°' 303 ] 

* [l-a(G)RA(l)TD(6, 10) 2 ]" 1 (3- 2 2d ) 

A ( 6 , 10)= FA64- - FA104- (3-22e) 

Here the variables are defined similarly to those of the 
preceding section. (Note that A(2, 4) 1 and A (2, 4 ) " have 

been combined in Equation (3-22) to give A(2, 4) and the 

definition of TD(6, 10) is similar in principle to that 

of TD(6, 8)). The absorbed insolation at the ground is 

given by 

I A1 0 ( 1 , 0) = FA10+ (1 - a (G) ) (3-23) 

(see Figure 5 ) . 

4 . Disposition of F(A) Insolation with a Low Overcast 
The set of equations which describe the behavior 
of the insolation with a low overcast deck are given 
below (see Figure 6) . 
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Figure 5, Schematic representation of F (A) insolation dis- 
position with an upper overcast layer only. 
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Figure 6. Schematic representation of F ( A ) insolation dis- 
position with a lower overcast layer only. 



35 



A(2, 


6)' = .27 IF (A) [U (2 , 6 ) Sec Z]° ,3 ° 3 


( 3 - 2 4 a ) 


A(6, 


8)' = .27 IF (A) [U (2 , 8 ) Sec Z)°* 3 ° 3 






n ^ n ^ 

- (U ( 2 , 6 ) Sec Z) * ] 


( 3-24b) 


FA8 + 


0 "3 0 1 

= F (A) [1-. 271 (U (2, 8 ) S ec Z) ’ ] 


( 3-24c ) 


FA9i 


= FA8+ [1-RA (2) -A(2) ] 


(3-24d) 


FA8t 


= FA8 4-RA ( 2 ) 


( 3 - 2 4 e ) 


A(6, 


8)" = FA8i [ .271 (U (6, 8)Sec Z)°* 3 ° 3 ] 


( 3 - 2 4 f ) 


A ( 6 , 


8) = A (6, 8) * + A (6, 8) " 


( 3 — 2 4g ) 


A ( 2 , 


o 3 n 3 

6)" = . 271FA8t [ (U (2, 8 ) Sec Z) 






— (U ( 6 , 8 ) Sec z ) 0 - 3 03 ] 


( 3- 24h ) 


CM 

< 


6) = A ( 2 , 6)' + A ( 2 , 6)" 


( 3 - 24 i ) 


A ( 8 , 


9) = FA8I A ( 2 ) 


( 3 - 2 4 j ) 
2 - 1 


FA10I 


= FA9 iTD ( 9 , 10) [1-a (G) RA (2) (TD (9, 


10 )) Z ] 

(3-24k) 


A ( 9 , 


10)= FA94- - FA10 + 


(3- 24£ ) 



These variables also have similar definitions to those in 
Section III-C-2. As before the absorption by the earth's 
surface is given by 

I Al 0(0, 1) = FAlOi (1 - a (G) ) (3-25) 



5 . Weighted F (A) Surface Insolation 

As previously described in the section on the 
F ( S ) insolation, the values of IA10 calculated for the 
cases (0, 0), (1, 1), (1, 0) and (0, 1) must be appro- 

priately weighted by the values described in Equation (3-12) 
in order to arrive at the composite F (A) surface-absorbed 
insolation. 
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IA10 = I AlO ( 0 , 0)W(0, 0)+ I A1 0 ( 1 , 1)W(1, 1) 



+ IAIO ( 1 , 0)W(l f 0 ) + I AlO ( 0 , 1)W(0, 1) (3-26) 

6 . Atmospheric Layer Absorption of F (A) 

In each of the four cases discussed, the F(A) in- 
solation was absorbed as it passed through the key atmos- 
pheric layers (2 f 6) and (6, 10). In the clear-sky case 

these absorptions are given by Equations (3-15) and 
(3-17). However, in the calculations of the cloudy-sky 
cases the various contributions have mer'ely been added 
together to give the total insolation absorbed in the 
layer (2, 6) and then in the layer (6, 10), even though 

cloudiness affects the radiation physics differently within 
the cloud and subcloud layers as Figures 3, 4, 5, and 6 

indicate . 

a. Overcast in Both Cloud Layers 

In this case the insolation absorbed in the 
two layers was calculated with values from the appropriate 
equations from (3-19a, ..., 3-19&) 



A(2 f 


6, i.i = 


A ( 2 , 


4) ’ 


+ A ( 2 , 


4) " 


+ A ( 4 , 


6) 


A (6, 




A ( 6 , 


8) 


+ A ( 8 , 


9) 


+ A ( 9 , 


10) 



where the subscript (1, 1) on A(6, 10) indicates the cloud 

layer combinations which will subsequently be weighted, 
b. Upper Overcast Only 

For the high level overcast only the insola- 
tion absorbed was calculated using the various equations 
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(3-22a, ...» 3 - 2 2 e ) . 



A(2, 6) = A(2, 4) + A(4, 6) 

A ( 6 , 10) = FA6 i - FA104- 

1 # u 



(3-28) 



c. Lower Overcast Only 

In the low level overcast cas<e the value of 
the insolation absorbed in the layers ( 2 „ 6) and (6, 10) 

was calculated using equations of (3-24a*..., 3-24£) 



A (2, 


6) 0.1 


= A ( 2 , 


6) ' 


+ A ( 2 , 


6) " 


(3- 


-29) 


A (6, 


10) 0,1 


= A ( 6 , 


8) ' 


+ A ( 6 , 


8 ) M + A ( 8 , 9 ) 


+ A ( 9 , 


. 10) 



d. Composite Weighted Absorption by Layers 

In this case, as in the previous subsection 
on weighted F ( A ) surface insolation, the values computed 
in each of the four cases must be appropriately weighted 
by the functions described in Equation (3-12) in order to 
arrive at a composite value of the F(A) insolation absorbed 
in the atmospheric layers (2, 6) and (6, 10). 



A ( 2 , 


6 ) =A ( 2 , 


6) o,o 


W (0, 


0) +A (2 , 


6) 1 1 W(l, 1) 






+ A(2, 


6) 1.0 


w (1 , 


0) +A (2 , 


6) o,i w(0 ' 15 


(3-30) 


A ( 6 , 


10) =A (6, 


10) o 


,o w ‘° 


, 0 ) + A (6 


, io) mi, i) 






+ A ( 6 , 


10) 1. 


0 WU, 


0) +A (6 , 


10 0,1 W(0 ' 15 


(3-31) 



e c The Absorptivity (ABA) by Layers 

In this study it was convenient to deal with 
the fractional absorptivity, or simply absorptivity, as 
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well as the actual insolation energies absorbed in each 
layer, given by Equations (3-32 and 3-33) and by Equation 
(3-26) at the ground. In the calculation of the absorp- 
tivity, the total incoming insolation at the top of the 
atmosphere must be known. This value was calculated us- 
ing Equation (3-32) 

FAD J = 2 . 00 ( R/R )~ 2 Cos Z (3-32) 

m 

Table II lists the instantaneous values of FADJ at .each 
gridpoint along the three data lines. The calculation of 
the f u 1 1- tr opospher ic absorptivity then follows from 



ABA 



A ( 2 , 6) 4 - A(6,- 10) 

^ FADJ 



(3-33) 



which is the ratio of the insolation absorbed in the tropo- 
sphere to the extraterrestrial insolation at the same time 
for each gridpoint. 

D. REFLECTION OF F (A) AND F(S) ENERGY 

In the calculations of the albedo and the fractional 
absorptivity of the ground both the F (A) and F(S) frac- 
tions of the solar insolation had to be considered. 

1 . Albedo 

That portion of the F(A) insolation which was re- 
flected back to space is given by 

REFA = [F (A) - A ( 2, 6) - A(6, 10) - IA10] (3-34) 
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TABLE II. Total extraterrestrial insolation ( FAD J ) at 
each grid-point of the three data lines (ly min”l). 



Line 1 


Line 2 


Line 3 


.899 


1.899 


1.346 


.923 


1.912 


1.374 


.947 


1.923 


1.403 


. 970 


1.933 


1.430 


. 995 


1.939 


1.458 


1.018 


1.946 


1.485 


1.042 


1.947 


1.510 


1.064 


1.947 


1.534 


1.085 


1.942 


1.558 


1.106 


1.934 


1.579 


1.125 


1.923 


1.597 


1.142 


1.907 


1.613 


1.158 


1.888 


1.626 


1.171 


1 „ 862 


1.634 


1.180 


1.834 


1.638 


1.187 


1.798 


1.637 


1 .190 


1.757 


1.629 


1.187 


1.710 


1.616 


1.181 


1.658 


1.593 


1.168 


1.598 


1.565 


1.149 


1.534 


1.526 


1.122 


1.462 


1.478 


1.088 


1.384 


1.419 


1.048 


1.301 


1.350 


.998 


1.210 


1.270 
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where the calculation equations of A(2, G> ) , A(6, 10), and 

IA10 have been described in Equations (3—30), (3-31), and 

(3-26) , respectively. 

The F ( S ) insolation reflected was simpler to calcu- 
late because the F(S) wavelengths do not undergo atmospheric 
absorption in the passage through the layers. The value 
of the reflected F(S) insolation is given by 

REFS = F ( S ) - IS10 (3-35) 

where IS10 was calculated by Equation (3—13). 

These values of REFA and REFS were then combined 
to give a value for REF, which is the total reflected energy 
passing outward through the tropopause (k=2) 

REF = REFS + REFA (3-36) 



The albedo is defined as the ratio of the total re- 
flected energy to the total incoming energy. This ratio 
was calculated by 



ALB 



REF 

FADJ 



(3-37) 



where FADJ is defined in Equation (3-32),. The values of 
REFA and REFS were computed at each grid-point. The re- 
sulting REF values are shown in cross-section form in 
Figures 11(a) and 11(b). Also, values of the albedo were 
computed. The global mean albedo over tlie 75 soundings 
turned out to be .412, based upon the cloud-cover para- 
meterization for C L ( 1 ) , CL ( 2 ) . 



4 1 



2 



Composite Absorptivity ( ABG ) by Earth 



The total insolation absorbed 
the sum of the weighted values of the 
portions of the incoming solar insola 

Q (E ) = I A10 + IS10 

Here Q( E) is the insolation absorbed 
face at each grid-point sounding. 

Using Q (E) , a fractional abso 
earth's surface was calculated. The 
the calculation was 



at the earth was 
F(S) and the F (A) 
ti am . 

(3-38) 

by tthe earth's sur- 

rptivity for the 
equation used for 



where the 
discussed 



ABG 



Q(E) 
FAD J 



calculations of Q(E) 



(3-39) 

and FAD J liave been previously 



3 • Check on Computations 

The actual computational process used in this model 
was checked against the condition that tlhe summation of the 
fractional values ALB , ABA , and ABG at each grid-point 
equals .96. This must be expected since we allowed for 
attenuation of the solar insolation in tlhe amount of four 
percent as it passed through the stratosphere. This check 
computation was satisfied at each grid-point. 



E. STATISTICAL COMPUTATIONS 
1 . Clear Sky Cases 

In order to verify the preceding calculations, 
statistical computations of several key values were obtained 
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using linear regression techniques from the BMD set of 
statistical programs (Dixon, 1973). The predictand of 
the first computation was ALB (0,0, Z), which is the 
albedo in the clear sky case. The predictor was the se- 
cant of the zenith angle (Sec Z). The best-fit equation 
which resulted was 

ALB ( 0 , 0 , Z ) = (-0.02726 + 0.11387 Sec Z) (3-40) 

with the regression coefficient R = .979. Equation (3-40) 
summarizes how the radiative model linearizes the effect of 
Sec Z in the clear sky case. The high R-value confirms 
the known result that Sec Z is a good predictor of the 
albedo (with clear skies), which is to be expected over 
an oceanic area. The mean albedo in the clear sky case, 
which corresponds to the average secant of the zenith 
angle, is as follows 

ALB ( 0 ~, 5T = .135 Sec Z = 1.428 

Two other regressions that were performed using 
the linear regression technique were ABA (0,0, M) and 
ABG (0,0, M) against the water-vapor mass path length M, 
where M is' given by 

1/2 1/2 

M = (U Sec Z) 7 [log (U Sec Z) 1 

This parameterization of the water vapor mass is essential- 
ly that of Hanson (1971), who used a similar absorber-mass 



43 



predictor in his empirical specification of ABA for both 
clear and partly cloudy skies. The resulting equations 
obtained with clear skies in this study for ABA (0,0, M) 
and ABG (0,0, M) were 

ABA (0,0, M) = .096 + ,031 (U Sec Z) 1</2 [lo^ 1Q (U Sec Z) 1 ^ 2 ] 

R = .984 (3-41) 

1/2 1/2 

ABG (0,0, M) = .800 - .035 (U Sec Z) / [logr 10 (U Sec Z) / ] 

R = .956 (3-42) 

The means and standard deviations of ABA (/0 , 0 , M) , ABG(0,0,M), 
and M, based on all of the original 75 soundings , are 
given below 

ABA (0, 0,M) = .137 O = .029 

ABG (070,M) = .687 O = .034 

-2 1/2 

M = 1.342 ( gm cm ) " a = .934 . 

2 . Relationship between Albedo in th:<e Cloudy and 
Clear Sky Cases 

As discussed previously, it was possible to compute 
the value of ALB in the clear sky case amd also in the com- 
posite cloudy sky case. In order to studi\y the relationship 
between these two cases a ratio ALB [ CL ( 1 ) CL ( 2 ) ] /ALB ( 0 , 0 ) 

for the same zenith angle at grid point (’I , J ) was formed. 

This value was used as the predictand ancl the total opaque 
c loud cover (CL) ? specified by the model was used as the 
predictor. Here CL is calculated using Equation (3-43), 
which approximates the effective cloud-cover by the two 
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layers of clouds: 



CL = CL ( 1 ) + CL ( 2 ) - CL ( X ) ( CL ( 2 ) (3-43) 

The regression equation which resulted was 

ALB [CL (1) ,CL (2) ] = ALB(0,0)(1 + 3.931CL - . 8 7 8 2CL 2 ) ( 3 -4 4 ) 

R = 0.956 



The 


means and standard 


deviations 


o f the 


values 


t ion 


(3-44) are 










ALB [CL (1) ,CL (2) ] 


= .412 


0 = 


..431 




ALB (0,0) 


= .135 


a = 


. 141 




CL 


= .680 


o = 


.742 



Note that the mean cloudy sky albedo generated by 
the model is 3.211 times that of the mean clear sky albedo. 
The high value of the correlation coefficient (R = .956) 
indicates reliability of the indicated dependence upon CL. 
This global mean albedo of .412 is somewhat higher than 
that recently reported [ALB = .300] by Raschke et al (1973). 
However, the mean value reported here does not have an 
input from the polar latitudes. Raschke, et al suggest 
that perhaps the reduced cloud amounts (CL) recently re- 
ported in low latitudes would account for the reduction 
compared to that computed here by the parameterization 
Equation (2-13), after Smagorinsky. 
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3 • Statistical Relationship between Ground Insolation 
in the Cloudy and Clear Sky Cases 

In this section as in the preceding one a ratio 
Y = ABG (CL) /ABG (0 , 0) was formed for each sounding case 

to determine the best-fit regression formula using CL as 
the independent variable. The formula which resulted was 

Y E ABG (CL) /ABG (0 , 0) = (1 - .66CL), R = .999 (3-45) 

The ratio of ABG ( CL ) /ABG ( 0 , 0) is essentially the 

corresponding ratio of the transmissions received by a 
pyrheliome ter . Quinn (1970) has summarized the empirical 
solar transmissivity formulas by many investigators and 
shows that some of the most reliable estimates are very 
nearly as formulated by Equation (3-45). For example, 
Budyko (1956) for the tropics lists 

ABG ( CL ) / ABG ( 0 , 0) = 1 - .655CL 

and Kimball (1928) lists the corresponding ratio as 

ABG (CL) /ABG (0 , 0) = 1 - .71C 

The mean statistics over the sample soundings gave 

ABG (CL) = .379 
ABG ( 0 , 0) = .688 

CL = .680 
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4 . Absorptivi t ies in Air 

Here the computed values of the ratio ABA (CL ) / 

ABA ( 0 , 0) were investigated by least squares with the 

result 

ABA (CL) /ABA (0, 0) = 1 + .37863CL, R = . 781 (3-46) 

with mean statistics for the 75-case population 

ABA "(CL)" = .169 

ABA (0 , 0) = .137 

CL = .680 

This indicates that the atmospheric absorptivity is some- 
what enhanced with the cloud-model solar dispositions in- 
troduced. An increasing solar absorptivity with increasing 
cloud cover is in general agreement with Plante's results 
(1972) in a one-layer cloud model. 

5 . General Conclusions 

The statistical relationships developed from the 
radiation model show consistency with results reported by 
other researchers. For example, Equation (3-41) gives 
ABA which closely fits Hanson's (1971) empirical results, 
and Equation (3-45) gives a ratio of transmittance which 
closely agrees with the empirical results of Budyko (1956) 
and Kimball (1928). Finally, the mean value of ABA in 
clear skies (ABA =.137) corresponds closely to the tropo- 
spheric mean absorptivity of 15 percent reported by London 
(1957). 
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IV 



TERRESTRIAL RADIATION 



A. THEORETICAL AND EMPIRICAL BASIS 

Sasamori (1968) developed empirical formulas for 
computing emissivities for long-wave flux calculations 
associated with the NCAR General Circulation Model* 

These empirical emissivity formulas are matched to the 
values developed by Yamamoto (1952) from radiative transfer 
theory, which have also been built into bis well-known 
Radiation Chart, This chart has proved to be quite accu- 
rate for computational purposes, and is used here as a 
schematic guide for integration of the radiative transfers 
through the various layers of interest in the soundings 
discussed in Section II e 

For utilization in the FNWC heating package the neces- 
sary emissivity formulas of Sasamori are those which will 
lead to expressions for 



net 


IR 


flux 


at 


earth 




net 


IR 


flux 


at 


level 


k = 6 


net 


IR 


flux 


at 


level 


k = 2 



* * 



* * 



The differences F^ - F, ^ and F^ - F ^ will then give the IR 
flux-divergence (i,e,, cooling rates) in the layers (6, 10) 

and (2, 6), respectively, A workable scheme for making 

such computations is presented below for various combinations 
of cloud cover [CL(1), CL(2)], 
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(4-1) 



For example, in order to compute 

B, 



F 10 = B 10 



7 



B = 0 



10 



£ (U, C, T) dB 

wc 



one must have a good representation of the emissivity £ 

wc 

as a function of both water vapor and CC> 2 absorber masses 

along the sounding, that is along the path of integration. 

Figure 7(a) illustrates the concept involved at level 

k = 10, with clear skies, where the hatched area represents 

the downward flux at k = 10. Here the radiative sounding 

is converted into emissivities as indicated in the figure*, 

The hatched area denoted by „ 

B 10 

£ dB 
wc 

B = 0 




is then the key quantity to be evaluated. For this integral 
evaluation, Sasamori (1968) proposes emissivity formula- 
tions of the type 



e wc ( 8, 10) = ) - 240 log 10 [U(Q, 10)] + .622 J 

+ .07262 |(1. - .62556 [U (8 , 10) + .0286]* 26 ) 

* [log 1() C ( 8 , 10) + 1.064] J (4-2) 



for temperatures T 210°K, where e is very nearly tempera- 
ture independent. The quantity in the first brace of Equa- 
tion (4-2) indicates the water vapor emissivity alone, 
whereas that in the second brace is the added contribution 
by the part of the CO^ emissivity transmitted through the 
water vapor "overlap" near 15 |im. 
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Equation (4-2) was used in the calculation of the atmos- 
pheric emissivities over any limits (6, ILtO ) , (4, 10), ... , 

(1, 10) for which running mass values U(6 ; , 10), C(6, 10), 

etc. were available from the radiative soundings. Hence 
Equation (4-2) is extremely versatile in its applicability. 

In order to complete the integration frorm level k = 1 to 
the origin (B = 0, £ = 1.0), the U(T^), C‘(T^) soundings 
were taken as isothermal at T = over the layer (0, 1) . 

However, it is still necessary, as Sasamori points out, to 
determine the hatched areas on the Yamamoto diagrams for 
the temperature-dependent emissivity £ wc E‘U(0, 10), T^]. 

This was computed after Sasamori as 

1 [U(0, 10), T ] = 8.34T * 353 1Og i0 U ( 0 ' 10) “ * 44 

wc 1 1 

*U(0, 10 ) 03455 ‘* 705 * [ 8 . 0/ ( . 3 5 3 log U(0, 10) + 3.56)] 

+ ,07262| (1. - .62556 [U(0, 10) + . 02186 ] * 26 ) 

*[log 10 C ( 0 , 10) + 1.064] J (4-3) 

The final term in the braces of Equation (4-3) involves 
essentially the same form of the CO^ correction to the 
total emissivity as appeared in Equation (4-2). The first 
term in Equation (4-3) is, however, temperature-dependent, 
as well as dependent upon the summed absorber masses rela- 
tive to the reference level (taken here as k = 10). It is 
clear that Equation (4-3) is applicable with appropriate 
final parameters U(0, k), C(0, k), and T , to the 
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determination of the curvilinear portion of the downward 
flux contribution through level k for all of the reference 
level cases to be considered below, e.g. # k = 2, 6. 

B. NET FLUX F* Q , WITH CLOUDY SKIES 

The radiative soundings have been evaluated in the 

form of the parameters U(k, 10), C(k, 10), T , level for 

K 

level. In addition there is listed for each radiative 
sounding at every grid-point the cloud-cover parameters 
CL ( 1 ) and CL(2), which in general are not both zero. The 
grid area may then be visualized as comprising the areal 
fractions W(0, 0), W(l, 1), W(l, 0), and W(0, 1) defined 

by Equation (3-12), of the following cases: 

(a) totally clear skies, (case 0, 0) 

(b) overcast in both layers, (case 1, 1) 

(c) overcast in upper layer only, (case 1, 0) 

(d) overcast in lower layer only, (case 0, 1) 

The resultant radiative effect of the true cloud cover 

over the grid area is then reconstructed by multiplying 

the results of the respective cases by the appropriate 

weight factors. It is expedient to discuss first the 

totally clear sky case (Figure 7(a)). 

* 

lo Net Flux F with Clear Skies 

If the hatched area F + ^^ °f Figure 7(a) is evalua- 
ted using the trapezoidal summation rule, then the net 

* 

flux F becomes 
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B(T) 



Figure 7(a). 
(case (0,0)). 



* 

Terrestrial net flux with clear skies, 

The unhatched area depicts 




B 
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(1-e ) 



wc 



B= 0 



dB 
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F 10 = B 10 - • 5{£ wc (8 ' 10)( V B 8 )+( V B 6 )[e W c (6 ' 10) 

+£ wc (8 ' 10)I + [e wc (4 ' 10)+£ wc (6 ' 10)1 [e wc (2 ' 10) 

+e (4,10) ] + [e (i,io)+e ( 2 , 10 ) ] ( b _ -b. ) + e [ (0,10) ,t, ] b } 
wc wc wc 2 1 wc l 1 

B k = 1.170403 * 10~ 7 T k ly(day) -1 (4-4) 

Here is the S t ephen-Bol tzman blackbody flux at T^ . 

* 

2. Net Flux F with Overcast Clouds in Both Layers 
Here Figure 7(b) is converted into a numerical 

form using the Sasamori type emissivities [Equations (4-2) , 
(4-3)] in the quadrature process. The result is 

F 10 U ' 1} = (B 10 - V [1 - - * 5£ wc (9 ' 10)] (4_5) 

where B is defined in Equation (4-4). 

* 

3. Net Flux F with Overcast Clouds in Upper Layer 

Only 1 

Here Figure 7 (c) is converted into a numerical 
form by the quadrature process after computing using 

Equations (4-2) and (4-3). The result is 

F io (l - 01 ■ <B 10-V (1 - - * 5te „ c (s ' 10 >< B 10-V 

+ [£„ (8, 10) +e (6, 10 ) ] (B -B ) ] } (4-6) 

wc wc 8 6 

The rationale of Equation (4-6) is that the downward flux 

is essentially B plus an increment which is in the upper 

6 
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12 4 6 8 10 

B(T> — *- 



Figure 7(b). Terrestrial net flux F with both high 
and low level overcast cloud decks, (case (1,1)). The 
unhatched area depicts 



* 



F 



10 




B 



10 
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wc 



) 
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1 2 4 6 8 10 

B(T) — <- 



* 

Figure 7(c). Terrestrial net flux F with high level 
overcast clouds in the layer k=4 to k=6 only, (case (1,0)) 
The unhatched area depicts 

B 10 

F* = / (1- e ) dB 

10 / wc 

B = 0 
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right corner of Figure 7(c), given by the summed combina- 
tion of the products of the form 



. 5 e * A B 
wc 



as these appear in the brackets in Equation (4-6). 

* 

4. Net Flux F^ q with Overcast in Lower Layer Only 
Here the numerical result follows immediately from 

the graphical depiction of Figure 7(d). This case becomes 
a replicate of the result of Section IV(B)2, with 

F* (0,1) = (B - B ) [ 1 . -,5e (9,10)] (4-7) 

10 10 9 wc 

* 

5. Composite F^ Q Calculations 

Here we use the principle enunciated in Section IV 

* 

(B) , and express F [CL(1), C L ( 2 ) ] as the properly cloud- 
weighted sum of the respective results of cases (a), (b) , 

(c), and (d). The result must have the form [see Equations 
(3-13 ) and (3-26) ] 



F 1q = W ( 0 , 0 ) 
+ W(1,0) 



10 (0 ' 0) 


+ W(l,l) 


P 10 U ' 1) 




* 




* 




10 (1 '°) 


+ W(0,1) 


H 

O 

O 


(4-8) 



with the cloud-weighted fractions, W(0,Q), W(l,l), W(1,0) 
and W ( 0 , 1 ) defined in Section III(B)3. The composite re- 
sult simplifies to 

P* 0 ■= d-CM2)l{ <B 10 — B 6 ) -.5 e uc (8,10) <B 10 - B 8 ) 

♦ + c wc (6.10)l (B a - B 6 )} + 
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B(T) 



Figure 7(d). Terrestrial net flux F with low level over- 
cast clouds in the layer k=8 to k=9 only, (case (0,1)). 

The unhatched area depicts 




B=0 



e ) dB 
wc 
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+ (1-CL(2) ) (l-CL(l) ) | B 6 -.5 [(e wc (6,10>+e wc (4,10) ) (B 6 -B 4 ) 

+ te wc (4,i0) + e wc (2,i0)] (B 4 -B 2 ) + [e wc (2,i t 0) + e wc (i,0)] (B 2 - Bl ) 
+ e wc C(0,10) ,T 1 )J J. +CL ( 2 ) | (b 10 -b 9 ) [!>.-. 5e wc (9, 10) ] J. 

(4-9) 

* 

Note that Equation (4-9) reduces to when CL(1) = 

* 

CL ( 2 ) = 0, and to F (1,1) when CL(1) = C.L ( 2 ) = 1.0. Simi- 

* 

larly the other two subcases of F^q may b)e recovered by 
setting up the weights corresponding to CL (1) = 1.0*, 

CL(2) = 0.0 and CL(1) = 0.0, CL(2) = 1 . 0 „ respectively. 

This procedure is similarly extended to levels k=6 
and k=2. The results are as listed in tine following sub- 
sections . 

* 

C. NET FLUX F WITH CLOUDY SKIES 
6 

Formulation of expressions for the urrhatched areas of 
Figures 8(a), (b), (c) and (d), with appropriate weight 

factors W(0,0), W(l,l), W(1,0), and W(0,1) applied, when 
summed lead at once to the following result: 

F* - (l-CL(D) {B e -.5t€ uc (6,8)(B 8 -B 6 ) +euc (4,6) (B 6 -B 4 ) 

* tC wc (4 '6) +E „c <2 ' 6)) (B 4 -B 2 )+re wc (2.6)+C wc (l,6)) (B 2 - Bl ) 

+ e wc ((0,6),T 1 )B 1 l| + (l-CL(l) ) (1-CL(2) ) j (B^-Bg) [1- 

. 5 (e (6 , 8 ) +e (6, 10)) ] 1 + CL (1) /(B -B,) [l-.5(£ (6,8))]) 

WC WC ) (86 WC J 

+ CL(1) (1-CL(2) ) <(B -B.) ( 1- . 5 [ £ (6,8) + £ (6,10) ] ) i 

( 10 8 wc wc ) 

(4-10) 
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Figure 8(a). 
(case (0 , 0) ) . 



* 

Terrestrial net flux F^, with clear skies 
The unhatched area depicts 



* 





( 1 - e ) dB 
wc 



B= 0 
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B(T) — 



* 

Figure 8(b). Terrestrial net flux F , with both high and 
low level overcast layers (case (1,1;). The unhatched 
area depicts 




B= 0 



e ) dB 
wc 



60 



1 2 4 6 ® 10 




* 

Figure 8(c). Terrestrial net flux F , with high level 
overcast clouds in the layer k=4 to R=6 only (case (1,0)). 
The unhatched area depicts 




B = 0 



(1 



e ) dB 
wc 
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1 2 4 6 8 10 




* 

Figure 8(d). Terrestrial net flux F , with low level 
overcast clouds in the layer k=8 to R=9 only (case (0,1)). 
The unhatched area depicts 




B = 0 



e ) dB 
wc 
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Equation (4-10) reduces to the properly formulated net 
fluxes for cloud cover cases (0,0), (1,1), (1,0), and (0,1) 

as these are depicted in Figures 8(a),(b),(c) and (d). 

* 

D * NET FLUX F 2 WITH CLOUDY SKIES 

* * 

Similar to F_. and F _ the net flux at level k=2 can be 
1 U 6 

* 

calculated using the trapezoidal rule. The net flux F^ 
is schematically depicted by the unhatched areas of Figures 
9(a) , (b) , (c) , and (d) . After these areas have been calcu- 
lated and weighted properly as described in the preceding 
two sections the following formula results 

F*= (l-CL(l) ) {B -.5 [e (2,4) (B,-B,)+(E (2 , 4) +e (2,6)) 

2 8 wc 4 2 wc wc 

<B 6- B 4 , + Ie wc <2 ' 6 > t£ „c (2 - 8)1 <V B 6» +e „c (1 ' 2) (B 2- B l' 

+e ( (0, 2) ,T. ) (B, )] ) +(1-CL(1) ) (1-CL12) ) 

WC 1 1 

{ (b — b_) [1-.5 (e (2,8)+e (2,10) ] } 

10 o wc wc 

+ CL ( 1 ) {b. - . 5 [ £ (2,4) (B -B,)+e (1,2)(B -B ) 

4 WC 4 2 WC 2 1 

+e ( (0,2) ,T n )B. ] } (4-11) 

WC 1 1 

This equation also reduces to the net fluxes for the four 
cloud cover cases (0,0), (1,1), (1,0), and (0,1) when the 

appropriate weight factors for the simplified cloud covers 
are set into Equation (4-11). 
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* 

Figure 9(a). Terrestrial net flux F , with clear skies 
(case (0,0)). The unhatched area depicts 




B = 0 
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wc 
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* 

Figure 9(b) 0 Terrestrial net flux F^, with both high and 
low level overcast cloud decks (case (1,1)). The unhatched 
area depicts 




B=0 



e ) dB 
wc 
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Figure 9(c). Terrestrial net flux with high level over- 
cast clouds in the layer k=4 to k=6 only (case (1,0)). 
The unhatched area depicts 
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wc 



6 6 



1 2 4 6 6 10 




B(T) 



* 

Figure 9(d) 0 Terrestrial net flux F , with low level 
overcast clouds in the lower layer, R=8 to k=9, only 
(case (0,1)). The unhatched area depicts 




B=0 



S ) dB 
wc 
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E. APPLICATION TO HEAT BALANCE COMPUTATIONS 

In computing the heat balance of the atmosphere, the 
layers (2,6) and (6,10) were considered to be undergoing 
solar heating (Section III), Now we require the long- 
wave cooling effects both instantaneously and averaged 
over the day (Section V). The long-wave cooling in (2,6) 
is symbolized by FU26 which is defined by 

* * 

FU2 6 = F^ - F _ 

2 6 

Likewise that in (6,10) is symbolized by FU610 where 

* * 

FU610 = F _ - F _ _ 

6 10 

is the flux-divergence in the lower half of the atmosphere. 

For computations of the radiative "balance" at the tropo- 
* 

pause only F of Equation (4-11) is required along with 

the net usable insolation. At the surface, a radiative 

balance is computed at each gridpoint involving the absorbed 

* 

insolation and the net flux F^ of Equation (4-9), the 
latter considered a negative contribution. As will be 
seen in Section V, the long wave flux divergences FU26 and 
FU610 both act as destabilizing influences in their respec- 
tive layers. 

F. STATISTICAL RESULTS 

1 . Downward Flux with Clear Skies 

Brunt (1932) demonstrated that the nighttime radia- 
tive measurements at the earth were expressible in a 
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statistical form 



F d • B 10 U + b/ * 1 



(4-12) 



where is the S tephen-Boltzman black body radiation 

4 

(B 10 = St*T 10 ) at the surface, e is the water-vapor 

pressure in mb, and a and b are the regression coeffi- 

* 

cients, and F. = B- _ - F__. The regression coefficients 
d 10 10 

referred to by Brunt were shown to be subject to wide 
variations depending on B and /e. 

In this study a statistical test of Equation (4-12) 
was carried out by linear regression using a program from 

the BMD set of statistical programs (Dixon, 1973) 0 The 

* 

predictand selected for these computations was F^, which 
included the downward contribution of CO^ as well as that 
of water vapor on the calculation of the net flux at level 
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Th 


* 




F „ ^ 


= B 


10 
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B _ „ 
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10 



10 



(4-13) 



The correlation coefficient of Equation (4-13) shows that 

/ — * 

and B-^q /e are good predictors of and therefore 

of F^ . Using the values of a and b in Equation (4-13) and 



solving for F^ we obtain 



F d = B 10 



or 



( l-a ) - b/e 



F, = B-,. (.69817 + .01314 /e) 

d 10 



(4-14) 
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The high value of the correlation coefficient indi- 
cates that Equation (4-13) may be regarded as an oceanic 

version of the Brunt F .-equation. Note however that most 

a 

of the dependence of is with and that the variation 

* / — 

of with B^ Q /e i- s relatively slight. ’This is as to be 

expected over the ocean, where nearly constant relative 

humidities might well exist over large areas. The mean 

values of F and B, ^ for clear-sky cases are as listed be- 
d 1 0 



F , = 637.09 ly day 
d 

843.47 ly day ^ 

so that F, averages 75.5 percent of B, - These statistical 
d 10 

averages are based upon the 57 grid-point soundings in the 

Northern Hemisphere only (those in the Southern Hemisphere 

have been deleted in order to eliminate tropical bias). 

2 • Statistical Relationship between Cloudy Sky and 
Clear Sky Cases 

Based upon the radiative sounding parameters listed 

in Table 1(b) and Equation (4-9) it was possible to compute 
* * 

F^ q [ CL ( 1 ) , CL ( 2 ) ] and F 1Q (0,0) based upon otherwise identi- 
cal parameters. 

For the purpose of this statistical study the ratio 

★ * 

F [ CL ( 1 ) ,CL(2)]/F^q( 0,0) was formed as the predictand, 
while an approximation to the composite total opaque cloud 
cover at the grid-point, defined in Equation (3-43) as 

CL = C L ( 1 ) + CL ( 2) ( 1-CL ( 1) ) 
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for each cloudy sky case was used as the predictor. The 
simple regression formula based upon a best-fit to the 57 
Northern Hemisphere data soundings became 

F 10 (CL) = F io {0 ' 0){1 " -7575CL) (4-15) 

with a regression coefficient, R = .983* A regression 

2 

analysis using CL as the predictor resulted in only a 
very small improvement of R. The hemispheric mean values 
of the various parameters in Equation (4—15) were as 
follows 



* 

F (CL) = 112.48 ly day 



F * 0 (°,°)= 206.38 ly day 1 (4-16) 

CL = 0.6103 



Here the superior bar means the sample mean was taken 
over the 57 Northern Hemisphere grid-points. 

The mean results of (4-16) indicate that the sur- 
face net flux is decreased by the ratio (93.9/206.38) = .455 

with CL = 0.6103. This is to say that the downward flux 

is on the average 45.5 percent larger than with the 
clear-sky case for the cloud-parameterization of Equation 
(2-13) . 

* * 

3. F 2 (CL) Relative to F 2 (0,0) 

Here the statistical parameter sought is the pre- 

dictand 

F* (CL) 

Y = -f (4-17) 

F 2 (0,0) 
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where both numerator and denominator are computed directly 
from Equation (4-11) with (a) the known CL(1), CL ( 2 ) used 
in the numerator, and with (b) CL(1) = CEj( 2) = 0 in the 

denominator. The regression result was cast in the desired 
form 

Y = 1 - A^CL)- - A 2 (CL) 2 

with A^ , A^ to be determined. The resulting regression 
equation was 

Y = 1 - .0742 (CL) - , 14 5 0 (CL) 2 (4-18) 

with a multiple correlation coefficient R * 0.9431, and a 
standard error of estimate .0529. The mean values of 
F* (CL) , F* (0, 0) , were 

F 2 (CL) = 453.19 ly day 1 
F* Q (0, 0) = 517.70 ly day" 1 
CL = 0.6803 

The mean CL here was based upon the 75 sounding-sample 
which contained a larger proportion of tropical cases than 
did those of the preceding section. 

It is to be noted that the infrared loss to space 
has not been greatly impeded by the assumption of normal 
cloud cover since no thick cirrus clouds have been in- 
cluded in the model. With deep cirrus the "greenhouse 
effect" would be increased for gridpoint columns so affected. 
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For comparison with satellite long wave flux to 

. -I 

space, which is usually given in ly min , the foregoing 

* -1 

F ^ values should be divided by 1440 min day 

Comparison with satellite measures of long wave 

* 

flux Should be made with F^t rather than with F^, which 
includes a downward flux contribution from the strato- 
sphere. For this purpose, computations were made with 

* 

F ^ of Equation (4-11) but omitting terms involving both 
levels " 2 " to "1" and " 2 M to "O" , which are always asso- 
ciated with downward flux (see Figure 9). The comparison 
of the grid-point sounding mean F [ CL ( 1 ) , CL ( 2 ) ] was then 
made with the NIMBUS III global average long wave flux 
to space for the appropriate date (after Raschke et al, 
1973). The comparison is as follows 



F 2 [CL(1) , CL ( 2 ) J = 
long-wave flux to space= 



.407 ly min 

.352 ly min 1 (after Raschke) 



The comparison is quite good considering that the results 
of Raschke et al are Northern Hemisphere means for a two- 
week period in April 1969, whereas the results of this 
study, .407 ly min is essentially based upon a mid- 

Pacific cross-section on April 25th, which tends to be 
warmer than the globe as a whole. 
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V 



MERIDIONAL CROSS-SECTIONAL DEPICTION OF RESULTS 
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Between these extremes all three meridional cross-sections 



were used in the calculation of the composite case. Fig- 
ure 11(a) and (b) depicts the composite cross-section with 
Figure 10 providing the key for locating the desired values. 

In the cross-section to be depicted the insolation 
rates at level k=2, (QAVG) are given as averages over the 
24-hour day (25 April 1973). This was done by employing 
Equation (3-1), but replacing Cos Z by Cos Z where 

Cos Z = [h Sincf) Sin6 + Cos<f> Cos6 Sin h]/7T (5-1) 

and h = hour angle = Arc Cos (-Tanij) TanS) 

<{) = latitude 

6 = declination (Smithsonian Tables (List, 1963)) 

Here the double bar over Cos Z represents a time mean over 
a 24 hour period, that is for a time period compatible 
with the long-wave flux transfers. 

At each latitude a value of the instantaneous absorp- 
tivities, reflectivities or transmissivities was known for 
the time 0000GMT, and for the purpose of the cross-section 
presented here, these were assumed constant even though 
the insolation rates changed over the day. The IR trans- 
fer rates and their flux divergences were assumed to be 
constant both for the day and the composite grid-point as 
determined from their 0000GMT soundings. Thus it was 
possible to compute: 
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a) the earth-atmosphere radiative balance at level 
k=2 by the formula 



QTROP = QAVG - (REFF * F ) (5-2) 



b) the radiative balance in the layer ( 2 , 6) as given 

by 

NET 2 6 = ABS26 + FU26 (5-3) 

c) the radiative balance in the layer (6, 10) as given 

by 

NET610 = ABS610 + FU610 (5-4) 

d) the radiative balance at the earth’s surface as 
given by 



NETSFC = ABSFC + F 



10 



(5-5) 



In addition to these mean grid-point items the cloud covers 
CL ( 1 ) / CL ( 2 ) and (CL) in parentheses are depicted in the 
applicable layer after averaging to the grid-point. At the 
tropopause, and at the earth's surface a global average of 
radiative balance or net was obtained by cosine weighting 
as follows 



Q = 



]£ Q (4> i ) cos <j>. 



(5-6) 



Cos (p 



where Q($^) is the radiative net value under consideration 

at latitude 6.. Some of these values are contrasted with 
i 
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the other climatological results cited below 



This Study Other Data 

1) global balance at 

the tropopause .048 •061, Raschke 

2) global balance at 

the surface .166 .140, Budyko 

The outgoing long wave radiation in these computations was 
also in excellent agreement with that of Raschke et al 
(1973) as noted in Section IV, The global balance of total 
radiation at the tropopause is less than that reported by 
Raschke et al for the same general time of year. However, 
it was noted that this result primarily reflects the 
smaller value of planetary albedo (ALB = .300) reported 

from satellite observations. 



Another set of results for comparison with those calcu- 



lated here are due to 


Budyko (1956) . 


Budyko calculated 


the radiation balance 


at the surface 


for 


ten-degree lati- 


tude bands. In his study, these values 


were taken for the 


balance at mid-points 


of these bands 


and 


compared with the 


computed output from 


the present model. 


The following 


table shows this comparison. 












Budyko 




This Study 




(Longitude 180°) 


5 “N 


. 158 




. 203 


1 5 0 N 


. 205 




.209 


2 5 °M 


. 236 




. 190 


3 5 °N 


. 201 




. 148 


4 5°N 


.119 




. 103 


55° N 


. 144 




. 064 
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These 



All of the 


above values are 


in 


cal cm 


2 . -1 
min 


comparisons 


would have been 


in 


closer 


agreement 


results of 


this study had been 


in the 


form of an 


average as 


Budyko's results 


were . 





if the 
annual 
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Figure 11(a). Mean Meridional Cross section, tropical section. 
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Figure 11(b). Mean Meridional Cross-section, higher latitude section,, 



VI 



CONCLUSIONS 



Using three meridional lines of soundings over the 
Pacific Ocean a radiative model for inclusion in the 
FNWC numerical weather prediction scheme was evaluated. 
These three meridians of data- sound ings vere used to 
compute a space-time average of the radiattional para- 
meters for one day over the Pacific Ocean.. This space- 
time average agreed fairly well with values derived from 
other studies, including satellite measurements. There 
were some minor differences which could be ascribed to 
the fact that the data used in this experiment was a one 
day space-time average while the comparison data was an 
average of several weeks for the satellite data, and 
over years in the case of surface climatological data 0 

Differences could also be induced by the parameteri- 
zation of cloud amounts as well as by the treatment of 
cloud reflectivity and absorptivity for solar wavelengths 
In this model, only two layers of clouds were considered 
and the same empirical values of the reflectivities and 
absorptivi t ies were used for all latitudes, which could 
cause discrepancies when compared with satellite data. 

Although the results of this data treatment agree 
fairly well with climatological comparisons cited at the 
end of Section V, refinements may be possible in the 
present model. This model should be tested using data 
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at different seasons of the year and over varying 



locales. With a full year of sampling, refinements in 
the cloud parameterization, atmospheric absorption 
values and reflection treatments could be developed 
for accurate use in the FNWC predictive model. 
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